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We  report  a  study  of  the  valence  band  dispersion  of  twisted  bilayer  graphene  using  angle-resolved 
photoemission  spectroscopy  and  ab  initio  calculations.  We  observe  two  noninteracting  cones  near  the 
Dirac  crossing  energy  and  the  emergence  of  van  Hove  singularities  where  the  cones  overlap  for  large  twist 
angles  (>5°).  Besides  the  expected  interaction  between  the  Dirac  cones,  minigaps  appeared  at  the 
Brillouin  zone  boundaries  of  the  moire  superlattice  formed  by  the  misorientation  of  the  two  graphene 
layers.  We  attribute  the  emergence  of  these  minigaps  to  a  periodic  potential  induced  by  the  moire.  These 
anticrossing  features  point  to  coupling  between  the  two  graphene  sheets,  mediated  by  moire  periodic 
potentials. 
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Much  effort  has  been  directed  toward  using  graphene  in 
electronics  and  optoelectronics  to  exploit  its  high  electrical 
conductivity  and  unique  Dirac  fermion  quasiparticles 
[1].  With  continuing  progress  in  fabricating  large-area 
graphene  sheets,  [2,3]  one  can  now  transfer  one  or  a  few 
graphene  layers  onto  desired  substrates  [4]  or  construct 
hybrid  multilayer  structures  [5,6].  Such  transfer  processes 
unavoidably  introduce  azimuthal  misorientation,  or  twist. 
Many  growth  processes  also  result  in  twisted  multilayers 
[7-9].  Envisioning  applications  involving  more  than  one 
graphene  sheet  for  specific  properties  [10-12]  therefore 
makes  it  important  to  understand  the  electronic  properties 
of  “twisted  graphene”  [13]. 

A  key  issue  is  the  electronic  interaction  between  twisted 
graphene  layers.  Theoretical  approaches  have  shown  that, 
for  twisted  bilayer  graphene  (TBG),  interlayer  interaction 
occurs  at  discrete  locations  within  the  Brillouin  zone  (BZ) 
[14-18].  Depending  on  the  twist  angle,  one  can  expect 
Fermi  velocity  reductions  or  the  emergence  of  van  Hove 
singularities  (vHs).  Transport  measurements  imply  that 
TBG’s  charge  carriers  near  the  Dirac  crossing  energy 
(. Ed )  behave  as  if  in  an  isolated  graphene  sheet,  confirming 
theoretical  predictions  for  a  large  twist  angle  [19,20]. 
Scanning  tunneling  microscopy  and  Raman  spectroscopy 
support  the  notion  of  interlayer  interaction  through  the 
presence  of  vHs  [21-23]  and  a  moire  [24].  On  the  contrary, 
angle-resolved  photoemission  spectroscopy  (ARPES) 
investigations  of  a  similar  system,  twisted  multilayer 
graphene  [i.e.,  >two  layers,  typically  grown  on  the  carbon 
face  of  silicon  carbide  (SiC)],  provided  no  evidence  of 
interlayer  interaction  across  the  entire  BZ  [25-27],  despite 
formation  of  moire  [28].  So  far,  ARPES  has  provided  little 
information  regarding  the  TBG’s  interlayer  interaction 
[29].  Thus,  questions  remain  on  the  existence,  extent, 
and  origin  of  its  interlayer  interaction  of  the  twisted  gra¬ 
phene  system. 


We  present  a  comprehensive  picture  of  electronic  dis¬ 
persion  in  TBG,  the  simplest  twisted  graphene  system, 
based  on  ARPES  and  density  functional  theory  (DFT) 
calculations.  We  observed  a  band  topology  consisting  of 
two  noninteracting  Dirac  cones  near  ED,  and  vHs  and 
associated  minigaps  away  from  ED,  where  the  two  layers’ 
Dirac  cones  overlap.  Our  experimental  results  provide 
unambiguous  evidence  of  the  interlayer  interaction  in 
TBG.  What  is  more,  we  observed  additional  minigaps  at 
the  boundaries  of  the  superlattice  BZ  associated  with  the 
moire  that  evolves  as  two  graphene  lattices  are  rotated  with 
respect  to  one  another.  Our  results  show  that  a  moire 
superlattice  gives  rise  to  a  periodic  potential,  altering  the 
electronic  dispersion  across  the  entire  BZ  according  to  its 
long-range  periodicity  and  not  just  where  the  states  from 
two  layers  overlap.  These  observations  illustrate  how  elec¬ 
tronic  dispersion  is  modulated  by  the  moire,  a  structure 
ubiquitous  in  superimposed  two-dimensional  (2D)  lattices 
(e.g.,  hybrid  multilayer  structures  [5,6]). 

We  fabricated  TBG  samples  by  transferring  graphene 
monolayers  grown  on  copper  foils  via  chemical  vapor 
deposition  [2,3,30]  onto  single-crystalline  epitaxial  gra¬ 
phene  monolayers  grown  on  a  hydrogen-terminated  SiC 
(0001)  (Si  face)  [31,32]  following  Ref.  [33].  This  fabrica¬ 
tion  procedure  results  in  >100  /xm  domains  with  random 
rotational  orientation  between  two  graphene  lattices. 
Within  each  domain,  the  twist  angle  is  relatively  constant 
[34].  Such  samples  allow  a  systematic  ARPES  study  of 
electronic  dispersion  primarily  on  a  single  domain  with 
minimal  effect  from  the  underlying  substrate  [35],  The 
underlayer’s  Dirac  cone  is  fixed  in  momentum  space 
( k  space),  while  the  overlayer’s  rotates  about  the  T  point 
of  the  first  primitive  BZ,  depending  on  the  twist  angle  9 
[33].  ARPES  measurements  were  conducted  at  beam 
line  7.0  of  the  Advanced  Light  Source  [36]  by  using 
95  eV  photons,  a  spot  size  of  —50  X  100  /mm2,  and  sample 
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T  —  100  K.  Given  the  photon  spot  size,  morphological 
variations  at  the  micron  scale  [23]  are  averaged  out  in  the 
ARPES  measurement.  Overall  energy  resolution  was 
—60  meV. 

DFT  calculations  were  conducted  by  using  VASP  [37] 
with  the  Ceperley-Alder  local  density  functional  [38],  as 
parameterized  by  Perdew  and  Zunger  [39],  in  the  projector 
augmented  wave  approximation  [40].  DFT  inherently  de¬ 
scribes  any  interlayer  electron  hopping  and  interaction.  We 
used  a  400  eV  plane-wave  basis  cutoff.  Correspondingly, 
optimization  of  single-layer  graphene  yielded  a  C-C 
separation  of  1.41  A.  Following  Shallcross,  Sharma,  and 
Pankratov  [15],  we  constructed  a  table  of  commensurate- 
moire  cell  sizes,  which  revealed  that  a  11.64°  twist  angle 
corresponds  to  a  TBG  supercell  with  a  repeat  distance  of 
8.54  A  containing  292  carbon  atoms  (146  in  each  layer). 
This  cell  corresponds  to  Shallcross’s  parameters  p  =  3  and 
q  =  17.  The  electronic  band  structure  at  this  twist  angle 
was  computed  for  comparison  to  the  ARPES  data  of  nomi¬ 
nally  0  =  —11.6°.  We  first  obtained  a  self-consistent  TBG 
charge  density  corresponding  to  a  9  X  9  equally  spaced 
sample  of  the  2D  superlattice  BZ  that  included  the  zone 
center.  We  then  computed  energy  levels  using  that  density. 
With  a  bilayer  separation  of  3.4  A,  local  density  approxi¬ 
mation  forces  on  carbon  atoms  along  the  bilayer  normal 
were  <0.02  eV/A. 

In  the  plots  shown  here,  calculated  DFT  data  are  shifted 
to  align  Ed  with  the  Fermi  level  ( EF ),  expanded  by  13%  in 
energy  to  account  for  many-body  interactions  [41],  and 
adjusted  for  the  average  doping  of  the  sample  [42].  The 
doping  level  (50  meV)  was  estimated  so  that  the  DFT  data 
best  match  the  ARPES  data  along  two  high-symmetry 
directions  [Figs.  2(d)  and  2(e)].  Unequal  carrier  concen¬ 
trations  between  two  graphene  sheets  were  not  considered 
in  the  calculation. 

The  electronic  states  of  TBG  can  be  described  by  two 
primitive  BZs  having  a  twist  angle  between  them.  The  case 


of  6  =  1 1.64°  is  illustrated  in  Fig.  1(a).  The  red  hexagon  is 
the  primitive  BZ  for  the  underlayer  graphene;  the  blue  one 
corresponds  to  the  overlayer  graphene.  The  two  small  half 
circles  overlaid  on  Figs.  1(a)  and  1(b)  are  the  measured 
photoemission  intensity  from  two  Dirac  cones  of  TBG  with 
8  =  —1 1.6°  [43].  These  intensities  are  the  constant  energy 
contour  at  0.4  eV  below  EF,  which  intersects  the  overlayer 
and  underlayer  cones  centered  at  the  K  and  Kg  points, 
respectively  [44].  Complete  separation  of  the  photoemis¬ 
sion  intensity  evidences  that  there  is  no  sign  of  interaction 
between  the  cones  close  to  ED. 

The  overlayer  and  underlayer  cones  exhibit  slightly 
different  spectra  as  seen  in  Fig.  1(c),  which  shows  photo¬ 
emission  spectra  along  the  line  connecting  K  and  Ke 
points.  Owing  to  photoelectron  attenuation,  the  underlayer 
cone  (left  cone  at  the  K  point)  displays  slightly  lower 
intensity.  Moreover,  although  the  overlayer  and  the  under¬ 
layer  graphene  sheets  are  both  p-doped,  their  earner  con¬ 
centrations  appear  slightly  different  as  evidenced  by  their 
dispersions  near  ED.  That  is,  the  size  of  the  overlayer 
Fermi  surface  (or  the  opening  of  the  cone  at  EF)  is  slightly 
smaller  than  the  underlayer’s  with  Ed  at  about  EP  + 
0.15  eV  and  EF  +  0.2  eV  for  the  over-  and  underlayer, 
respectively.  This  small  shift  in  ED  is  attributable  to  the 
smaller  influence  of  the  substrate  on  the  overlayer,  a  result 
of  screening.  Correspondingly,  epitaxial  graphene  on 
hydrogen-terminated  SiC  is  slightly  p-type  [32,45], 

Unlike  near  ED,  interactions  between  the  two  Dirac 
cones  are  observed  at  higher  electron  binding  energy. 
Figures  2(a)-2(c)  show  photoemission  intensity  patterns 
at  EF  —  0.8,  — 1.0,  and  — 1.3  eV,  where  the  contours  of  the 
cones  deviate  significantly  from  the  expected  circular-to- 
triangular  band  topology  of  monolayer  graphene  [46].  The 
superlattice  BZ  (black  hexagon)  and  its  high-symmetry 
points  Ts,  Ks,  and  K's  are  overlaid  on  these  experimental 
curves  to  highlight  the  superlattice’s  influence  on  the  in¬ 
terlayer  interaction.  [Note,  throughout,  that  subscript  “s” 


FIG.  1  (color  online),  k  space  representation  of  TBG  with  9  =  —  1 1.6° .  (a)  Photoemission  intensity  contour  of  the  two  Dirac  cones  at 
the  electron  energy  EF  —  0.4  eV,  and  the  primitive  BZs  of  the  underlayer  and  overlayer  (red  hexagon  including  the  K  point  and  blue 
one  including  the  Kg  point).  Darker  shades  indicate  higher  photoemission  intensities.  The  small  black  hexagon  is  the  moire  superlattice 
BZ  of  the  (p,  q)  =  (3,  17)  [15]  commensurate  TBG.  (b)  Enlarged  image  of  (a)  near  the  two  cones,  (c)  Photoemission  spectra 
intersecting  two  cones  at  K  and  Kg  points.  The  red  (left)  and  blue  (right)  dashed  lines  illustrate  the  under-  and  overlayer  cones, 
respectively,  obtained  by  extrapolating  the  measured  dispersions  [62]. 
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FIG.  2  (color  online).  Electronic  dispersions  of  the  two  inter¬ 
acting  Dirac  cones  (9  =  ~11.6°).  (a)— (c)  Photoemission  inten¬ 
sity  contours  at  EF  —  0.8  eV  (a),  EF  —  1.0  eV  (b),  and 
EF  —  1.3  eV  (c).  Black  hexagons  (thick  line)  indicate  the  moire 
superlattice  BZ  of  the  commensurate  TBG.  T(,  Ks,  and  K's  are 
among  its  high-symmetry  points.  K  and  Kg  points  are  both  Ks 
points  in  the  superlattice  BZ.  Green  hexagons  (thin  line)  are 
minizones  of  a  continuum  model  with  a  Dirac  point  at  its  zone 
center,  (d)  Photoemission  spectra  and  the  DFT  states  bisecting 
the  two  cones  and  (e)  the  one  orthogonal  to  (d).  Their  directions 
are  indicated  in  (c)  by  horizontal  and  vertical  black  arrows.  The 
schematic  to  the  right  of  (e)  shows  the  orientations  of  the 
photoemission  patterns  relative  to  the  two  primitive  Dirac  cones 
without  interaction.  DFT  states  are  shown  as  (blue)  dots. 
Calculated  states  matching  the  ARPES  data  are  highlighted  by 
blue  (thick)  and  green  (thin)  circles. 

indicates  high-symmetry  points  associated  with  the 
superlattice  BZ.]  Note  that  while  the  cones  exhibited 
“monolayerlike”  topology  near  ED  [Figs.  1(a)  and  2(a)], 
at  higher  binding  energy,  the  two  bands  merge  near  the  K's 
point  [Fig.  2(b)].  This  is  a  first  indication  of  their  interac¬ 
tion.  At  this  intersection,  nested  parallel  bands  emerge  to 
the  left  of  the  cones  [i.e.,  towards  the  origin  in  k  space  as 
indicated  by  a  red  arrow  in  Figs.  2(b)  and  2(c)],  which 
exhibit  an  anticrossing  behavior.  This  same  behavior  is  not 
seen  towards  the  right  [cf.  the  blue  arrow  in  Fig.  2(b)] .  The 
reason  is  addressed  below. 

Emergence  of  the  anticrossing  of  the  two  bands,  or 
minigap  formation,  results  from  coupling  between  the 
two  Dirac  cones.  To  illustrate.  Fig.  2(d)  displays  the  pho¬ 
toemission  spectra  along  the  horizontal  black  arrow  in 
Fig.  2(c),  which  bisects  the  two  cones.  Note  that  the  tt 


state  is  split  around  the  K's  point  (cf.  the  red  arrow).  This 
splitting  is  also  seen  in  the  DFT  electronic  levels  [blue  dots 
in  Fig.  2(d)].  The  anticrossing  behavior  can  be  understood 
in  terms  of  vHs,  when  the  orthogonal  direction  [vertical 
black  arrow  in  Fig.  2(c)]  is  examined.  Figure  2(e)  shows 
photoemission  spectra  and  DFT  results  along  this  direc¬ 
tion,  where  the  upper  “A7” -shape  and  the  lower  inverted 
“V” -shape  bands  correspond  to  the  left  and  right  nested 
parallel  bands  in  Fig.  2(c),  respectively.  By  noting  that  the 
M-shaped  band  in  Fig.  2(e)  is  the  same  as  the  upper  split 
state  in  Fig.  2(d),  it  is  apparent  that  these  states  have  both 
positive  and  negative  masses,  creating  a  saddle  point.  Thus, 
as  a  consequence  of  coupling  between  the  two  layers’ 
cones,  vHs  occur  at  the  anticrossing. 

Besides  the  vHs,  faint  states  reside  within  the  minigap 
near  the  red  arrows  in  Figs.  2(c)  and  2(e);  however,  they  do 
not  appear  in  the  DFT  calculation.  We  postulate  that  they 
are  due  to  the  areas  where  the  interaction  between  two 
layers  is  reduced  within  a  TBG  domain.  Such  locations  are 
attributable  to  topographical  defects  like  ripples  and 
blisters  [47].  Low  energy  electron  and  atomic  force  micro¬ 
graphs  support  their  presence  on  a  length  scale  much 
smaller  than  our  photon  spot. 

The  photoemission  intensity  contours  shown  in  Fig.  2 
include  an  additional  interacting  feature  not  explained  by 
direct  interaction  of  the  two  layers’  Dirac  cones.  The  green 
arrows  in  Figs.  2(b)  and  2(c)  highlight  a  splitting  in  the 
overlayer  cone  around  the  K's  point,  along  a  direction 
extending  into  the  upper-left  superlattice  BZ.  For  more 
details,  we  take  a  second  derivative  of  the  photoemission 
intensity  with  respect  to  electron  energy,  as  shown  in  Fig.  3. 
Red  and  blue  circles  in  Fig.  3(a)  highlight  under-  and 
overlayer  cones  and  help  illustrate  that  the  new  feature 
appears  not  as  a  consequence  of  these  cones’  intersection 
but  because  of  the  presence  of  a  “new”  cone  centered  on 
the  moire  superlattice  K' s  point  (black  circle).  Its  disper¬ 
sion  is  displayed  in  Fig.  3(b),  along  a  line  from  this  new 
(black)  cone  to  the  overlayer  (blue)  cone  [i.e.,  the  green 
arrow  in  Fig.  3(a)].  Similar  to  the  vHs  observed  in  Fig.  2(e), 
an  additional  vHs  is  observed  in  both  the  ARPES  and  DFT 
results  in  Fig.  3(b).  We  attribute  this  new  cone  and  the 
additional  vHs  forming  along  with  it  to  adiabatic  umklapp 
scattering  in  the  superlattice  periodic  potential  [48,49]. 
They  could  not  be  present  if  the  electrons  of  one  layer 
were  not  responding  to  the  periodic  potential  imposed  by 
the  other,  thus  confirming  that  the  two  graphene  layers  are 
not  isolated  but  sense  each  other. 

The  ramifications  of  the  periodic  potential  applied  to 
graphene  (in  this  case  induced  by  a  moire  superlattice) 
should  have  intriguing  consequences  [50-52].  In  “normal” 
2D  materials,  applying  a  periodic  potential  results  in  the 
isotropic  opening  of  the  minigap  over  the  entire  boundary 
of  the  minizone  defined  by  the  potential’s  periodicity. 
Graphene’s  response  is  quite  different  because  of  the 
chiral  (pseudospin)  nature  of  the  wave  functions. 
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FIG.  3  (color  online).  Second  derivative  of  the  ARPES  inten¬ 
sity  with  respect  to  the  energy  (9  =  ~11.6°).  (a)  Contours  at 
EF  —  0.8  eV.  The  black  hexagons  are  the  superlattice  BZ.  The 
red,  blue,  and  black  circles  with  “U”,  “O”  and  “AT”  illustrate 
the  locations  of  underlayer,  overlayer,  and  moire  superlattice 
Dirac  cones,  respectively,  (b)  Processed  photoemission  pattern 
along  the  green  arrow  in  (a)  and  DFT  states  (blue  dots)  con¬ 
necting  K' S-K' S-K' ,  points.  Circles  (green)  highlight  DFT  states 
matching  the  ARPES  data. 

Accordingly,  the  periodic  potential  does  not  open  the 
minigap  along  the  entirety  of  the  minizone  boundary 
but  only  at  certain  locations.  Thus,  moving  along  the 
minizone  boundary,  gaps  will  emerge  and  disappear. 

To  examine  this  effect,  following  Park  et  al.  [51],  we 
define  the  minizone  [green  hexagons  in  Fig.  2(a)]  by  trans¬ 
lating  the  superlattice  BZ,  so  the  center  of  the  superlattice 
BZ  matches  K  and  K0  points.  The  line  connecting  K's  and 
T,  points  is  along  the  minizone  boundary  [53].  In  this  view, 
coexistence  of  band  splitting  and  crossing  [shown  by 
ARPES  and  DFT,  red  and  blue  arrows  in  Fig.  2(d)]  is  a 
consequence  of  the  periodic  potential  induced  by  the  moire 
superlattice.  Hence,  the  nonconstant  gap  occurs  because  of 
graphene’s  chiral  wave  functions. 

Supporting  this  conclusion,  the  periodic  potentials  of  a 
moire  should  vary  on  a  much  longer  length  scale  than  the 
interatomic  distance.  Thus,  a  slight  shift  of  one  graphene 
sheet  relative  to  another  should  only  affect  the  electronic 
dispersion  of  TBG  weakly.  DFT  calculations  involving 
translations  of  one  of  the  two  graphene  sheets  by  a  fraction 
of  interatomic  distance  confirm  this.  We  conclude  that 
TBG’s  electronic  dispersion  evolves  from  two  rotated  gra¬ 
phene  sheets  subject  to  a  long-range  potential  of  the  moire 
superlattice  evolving  between  them.  Alternatively,  TBG 
comprises  two  graphene  sheets,  each  subject  to  a  periodic 
potential.  This  provides  a  simple  way  to  understand  many 
of  the  unique  features  alluded  to  in  previous  theoretical 
studies  [51,52]. 

Incidentally,  the  additional  interacting  state  does  not 
appear  at  the  underlayer  cone  highlighted  by  the  red  circle 
in  the  data  presented  but  did  appear  in  other  data  for 
different  (typically  smaller)  twist  angles  [54].  If  there  is 
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FIG.  4  (color  online).  Photoemission  intensity  patterns  (a),(c) 
and  EDCs  (b),(d)  displaying  the  minigap  as  a  function  of  9.  (a), 
(b)  9  =  ~5.6°,  (c),(d)  9  =  ~12.0°.  The  photoemission  pattern 
bisects  the  two  cones  similarly  to  Fig.  2(d).  (b)  and  (d)  are  EDCs 
at  the  black  lines  in  (a)  and  (c),  respectively,  fitted  to  Voigt 
functions  [thin  (red)  lines]. 


yet  another  new  Dirac  cone  present  in  Fig.  3(a),  we  expect 

it  to  be  centered  on  kx  —  1.4  A-1,  ky - 0.34  A-1  and 

have  the  band  topology  similar  to  the  new  cone  highlighted 
by  the  black  circle.  Although  there  is  a  state  near  where  we 
expect  to  see  the  new  cone,  its  shape  is  quite  different.  We 
therefore  suspect  that  the  data  at  the  lower  ky  in  Fig.  3(a) 
originate  from  another  TBG  domain  having  a  slightly 
different  twist  angle. 

Regions  of  AB  stacking  in  the  moire  superlattice 
dominate  the  interlayer  interaction  for  twist  angles  >5° 
[18],  resulting  in  the  minigaps  seen  in  Fig.  4.  Figures  4(b) 
and  4(d)  show  the  energy  distribution  curves  (EDCs)  half¬ 
way  between  the  K  and  Kg  points.  For  twist  angles  of  —5° 
to  —12°,  the  energy  separation  (peak-to-peak)  stays  near 
—0.2  eV  (cf.  the  red  arrows).  This  value  is  of  the  same 
order  of  magnitude  as  the  interlayer  interaction  parameter 
of  Bernal  bilayer  graphene,  —0.4  eV  [55].  We  attribute  the 
relatively  unvarying  magnitude  of  the  minigap  with  twist 
angle  to  the  persistence  of  local  AB  stacking  within  the 
moire  [56].  The  large  real-space  moire  superlattice  ensures 
the  existence  of  AB  stacking  for  all  twist  angles. 

Last,  we  offer  plausible  rationales  for  the  absence  of  some 
of  our  DFT  energy  levels  in  the  corresponding  ARPES  data 
[see  small  blue  dots  in  Figs.  2(d),  2(e),  and  3(b)].  First,  the 
structure  factor  associated  with  ARPES  may  preferentially 
increase  the  intensity  of  certain  states.  This  has  been 
observed  in  measurements  wherein  intensities  are  strongly 
enhanced  when  a  surface  state  overlaps  a  bulk  state  in  k 
space  [57].  Following  this  argument,  we  presume  that  the 
TBG  states  overlapping  those  of  a  noninteracting  graphene 
sheet  would  appear  strongly  in  the  ARPES  measurement. 
Consequently,  the  measured  photoemission  intensity 
matches  only  a  small  subset  of  the  DFT  calculated  states. 
Disorder  in  the  TBG  including  mechanical  distortions 
provides  another  possibility.  As  we  saw  via  low  energy 
electron  diffraction,  the  twist  angle  in  our  samples  varied 
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slightly,  over  a  few  micrometer  length  scale  [33,34,58]. 
Because  of  the  small  superlattice  BZ,  the  experimentally 
observed  states  from  slightly  different  twist  angles  would 
be  broadened  with  their  intensity  decaying  rapidly,  espe¬ 
cially  for  those  created  by  folding  at  superlattice  BZ 
boundaries.  The  likelihood  of  observing  these  folded  states 
would  decrease  correspondingly. 

Coupling  between  the  electronic  states  and  the  superlat¬ 
tice  periodic  potential  have  important  implications  for 
twisted  multilayer  graphene  and  hybrid  2D  multilayer 
stacks.  Based  on  our  study  of  TBG,  the  superlattice  BZ 
of  a  multilayer  graphene  (>  three  layers)  is  expected  to  be 
smaller  (thus  longer  periodicity  in  real  space).  Previous 
theoretical  work  has  shown  that,  with  an  increase  in  spatial 
period,  the  apparent  minigap  shrinks  [51],  leading  to  ef¬ 
fectively  noninteracting  states.  This  is  consistent  with  re¬ 
ported  experimental  results  [27].  Second,  any  hybrid 
multilayers  based  on  transferring  2D  materials  will  unavoid¬ 
ably  induce  moire  superlattices  and  thus  subject  the  system 
to  a  periodic  potential.  This  potential  influences  the  disper¬ 
sion  and  thus  the  properties  of  the  multilayer  stack. 
Although  transfer  techniques  now  offer  the  possibility  of  a 
wider  class  of  2D  materials  [59,60]  much  as  heteroepitaxial 
growth  does  [61],  understanding  how  these  layers  change  as 
they  are  stacked  together  and  mutually  interact  is  prerequi¬ 
site  to  leveraging  their  properties. 
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